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Modern high-speed computers are capable of assisting in the 
reduction and storage of data obtained from gas -liquid partition 
chromatography (GLPC) measurements. In the past, in order that a 
computer could be used, it was necessary to prepare machine language 
programs. This difficulty can now be largely avoided by the application 
of FORTRAN programing, 1 since by this technique the computer is used 
to aid in the preparation of the program. To illustrate this approach, a 
typical set of experimental GLPC data for a homologous series of 
alcohols is reduced by a FORTRAN program developed at this Labora- 
tory. The program is designed to give values for the specific retention 
volumes, the partition coefficients, and the number of theoretical plates 
for each peak appearing on a chromatogram. In the course of the data 
reduction other parameters of importance to the interpretation of results 
from GLPC measurements are also obtained. 


iThis paper presents the results of one phase of research carried out at 
the Jet Propulsion Laboratory, California Institute of Technology, under 
Contract No. NAS 7-100, sponsored by the National Aeronautics and Space 
Administration. 
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The reduction of GLPC data by FORTRAN is a two step process. 
The first step requires the preparation of a series of FORTRAN state- 
ments that define the mathematical operations by which the GLPC data 
are to be analyzed; this constitutes the source program. A FORTRAN 
Compiler translates the source program by means of punched cards 
into the object program, which is a set of elementary instructions that 
the computer can understand. The object program is placed into the 
memory of the computer. The second step is to place the experimental 
measurements taken from the chromatogram on punched data cards, 
upon which the computer carries out the designated mathematical 
operations. The computer presents the reduced data on either IBM 

cards or as a typed output. Further details for the preparation of the 

2 

FORTRAN statements are adequately documented elsewhere. 

Table I gives the source program and symbolism used for the 
reduction of the GLPC data. This program is intended to be used with 
an IBM 162 0 Data Reduction System, a 1622 Reader Punch, and a 162 0 
FORTRAN Format Compiler.' Although, it can be modified for use with 
other equally suitable data reduction systems. The program represents 
in FORTRAN language the mathematical expressions recommended by 
the ASTM Committee D-2 on Petroleum Products and Lubricants for 
the calculation of the corrected retention volume (VpL the specific 


*IBM Corporation, 1271 Avenue of the Americas, New York 20, New York. 


retention volume (V g ) measured at column temperature and ambient 

pressures, and the number of theoretical plates. The program also 

contains expressions for the corrected retention volume V g of Littlewood; 4 

the specific retention volumes measured at either column temperature 

(V 1 760) or 0°C (V 760) but corrected to standard pressure; the partition 
g g 

coefficient K ? gQ, measured at column temperature and corrected to 
standard pressure; and the pressure correction factor (f) of James and 
Martin. ^ 
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Table I. Source Program and Symbolism for GLPC Data Reduction 


A. Source Program 


34.35 


1 FORMAT (15X.4QHGAS LIQUID CHROMATOGRAPHY RETENTION DATA) 

2 FORMAT (//l6x p 2HFU,12x.3HPIG p 11X, 1HP,13X,IHT»13X, IHZ) 

3 FORMAT (10x p 3HSWT # nx,iHRTAIR.9X.lHD # J3X t 2HTC» J2X* iHN) 

4 FORMAT (10X,4HFSTP i IOX t 2HPT i lix # 2HPI # 12X # 2HPR,l2X,2HPF) 

5 FORMAT (10X.4HPGCF.10X P 2HVP,12X,2HTR) 

fa FORMAT MOx 2HRT 12 x 1 HW 1 3X 4HVGP7 , 1 OX ,3HVR0 p 11 X , 3HVG7 ) 

7 FORMAT ( 10x:3HCRT, n x* 2HtP, I ix ,3HVGP, l 1 X.3HCAY , 1 l X , 2HVG/ ) 

8 FORMAT (9F6.0.I2) 

9 FORMAT (5F14.4) 

19 FORMAT ( 4f 14 .4 ,9X » 13 ) 

10 FORMAT ( 5F 1 4.3 ) 

20 FORMAT (3F14.3) , y A c , n . 

11 FORMAT (F6.0 p F6.Q ! F6.0,F6.Q,F6.0 p F6.Q,F6,0,F6.0/ 

12 FORMAT (5F14.3/5F14.3/) 

IF (SENSE SWITCH 3} 30,31 

30 PRINT 1 

31 PUNCH 1 

00 90 1-1,100 

READ 8,FU p PIG,P p T.Z,SWT,RTAIR p D p TC p N 

IF (SENSE SWITCH 3) 32 p 33 

32 PRINT 2 

33 PUNCH 2 
IF (SENSE SWITCH 3) 

34 PRINT 9 p FU p PIG p P p T, 

35 PUNCH 9 p FU p PIG p P p T Z 
IF (SENSE SWITCH 3) 36,37 

36 PRINT 3 

37 PUNCH 3 

IF (SENSE SWITCH 3) 64,65 

64 PRINT 19 P SWT P RTAIR P D,TC.N 

65 vP-EXP(2:3O2jb5*(»:i67^^17S0.28fc/(235.O + T)))) 

FSTP-(FU*(P-VP)*273. 1 )/( (T+273 .1 )*7oO.) 

PT- ( P*1 4 .696) /76G . 
p|®( (PI G+PT >*760. )/1 4.696 
PR*P I /P 
PF=PR**2-1. 

PGCF = 1 ,5*PF / ( PR**3-1 . ) 

TR-1.0E4/(TC+273.1) 

IF (SENSE SWITCH 3) 3b, 39 

38 PRINT 4 

39 PUNCH 4 . , . , 

IF (SENSE SWITCH 3) 42,43 

42 PRINT 10,FSTP,PT,PI ,PR,PF 

43 PUNCH 10,FSTP,PT,PI ,PR p PF 
IF (SENSE SWITCH 3) 40,41 

40 PRINT 5 

41 PUNCH 5 ,, 

IF (SENSE SWITCH 3) 66,67 

bfa PRINT 20 , PGCF , VP , TR 
67 PUNCH 20 , PGCF p VP p TR 

IF (SENSE SWITCH 3)44,45 

44 PRINT 6 

45 PUNCH 6 

IF (SENSE SWITCH 3)46,47 
4b PRINT 7 

47 PUNCH 7 

DIMENSION CM( 100) , X ( 100) 

K-1 , 

21 READ 1 1 ,CM(K) ,X(K) P CM(K+1 
K-K+4 

IF (K-N) 21,21,22 

22 DO 80 J-1 , N 
RTN-RTA I R/Z 
RT-CM(J)/Z 
CRT c RT-RTN 
W“X ( J ) / 7- 

VG“FPP*273 . 1 / (T+27 3.1 ) 

VG7-VG*P/7bO. . , 

VGP“FPP*(TC+27 3 .') /(t+27 3.1 ) 

VGP7“VGP*P/780. 

CAY-(VG7*(TC+273.I )*D J/273.1 
TP-16. *((RT/W)»*2) 

4b PR I HT EI 1 ^RT^wI VGP7 , VRO ! VG7 .CRT , TP . VGP , CAY . VG 
49 PUNCH 1 2 ,RT ,W, VGP7 , VRO , VG7 , CRT , TP , VGP , CAY , VG 
bO CONTINUE 
90 CONTINUE 
END 


,X(K+I ) ,CM(K+2) ,X(K+2) ,CH(K+3) ,X(K+3) 
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B. FORTRAN Source Program Measurements 


FU Flow measured with a soap bubble meter (ml per min). 

PIG Pressure inlet gauge (lb per square in.). 

P Barometric pressure (mm of mercury). 

T Temperature at which gas flow measured (°C). 

Z Chart units per min (may be centimeters or inches as long 

as consistent units for retention time and peak widths are used). 

SWT Weight of liquid phase in column (grams). 

RTAIR Retention time of air in same units as Z. 

D Density of liquid phase at column temperature (grams per ml). 

TC Column temperature (°C). 

N Number of peaks in the chromatogram which are computed. 

CM Retention distance on chart in units which are consistent with Z. 

X Peak width in units consistent with Z. 

C. Column Parameters from Reduction of Data 

FSTP Flow reduced to standard temperature and pressure (ml per min). 
PT Atmospheric pressure (lb per square in.). 

PI Inlet pressure (mm of mercury). 

PR Ratio of inlet to outlet pressures. 

PF Pressure factor. 

Pressure gradient correction factor. 
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PGCF 


VP Vapor pressure of water at temperature T (mm of mercury). 

RT Reciprocal of absolute column temperature times ten thousand. 

D. Reduced Data for Components on Chromatogram 

RT Retention time of component (min). 

W Peak width of component (min). 

CRT Corrected retention time (min). 

TP Theoretical plates. 

CAY Partition coefficient. 

VRO Corrected retention volume (ml). 

VG Specific retention volume at zero degrees centigrade but not to 
standard pressure (ml per gram). 

VG7 Specific retention volume at standard temperature and pressure 
(ml per gram). 

VGP Specific retention volume at column temperature and at ambient 
pressure (ml per gram). 

VGP 7 Specific retention volume at column temperature corrected to 
standard pressure (ml per gram). 
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The data cards should be punched with nine fields of six columns 
each with one field of two columns at the end. The column input measure- 
ments are read by the computer in the following order: FU, PIG, P, T, 

Z, SWT, RTAIR, D, TC, and N, The first nine fields are floating point 
numbers with N being a fixed point number. If N is less than ten it must 
have a zero to fill the field. The input measurements for the components 
are read as eight fields of six columns each with four pairs of retention 
distances and peak widths per data card. Any fields not used must be 
completed with floating point zeros. In some instances it may not be 
practical to measure peak width. In these cases, the best procedure is to 
repeat the value for the retention distance as the peak width and this leads 
to the number 16 which appears as a fictitious number of theoretical 
plates. 

Figure 1 shows a chromatogram for a series of homologously related 
alcohols, and presents the proper organization of the input measure- 
ments for placement onto a data card for setting up the program. In 
Table II the computed data for the homologous series of alcohols is given so 
that it may be used as a check case. The time required to compute the 
check chromatogram is about 35 seconds. 

To make the program economical of computer time there is a 
provision which allows the computed data to be punched on IBM cards. This 
is controlled by sense switch 3 on the 1620 console. By having sense 
switch 3 in the on position, both typewriter and punch card output of data 
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results. This mode is slow, but allows checking at the start of the run. 
When sense switch 3 is in the off position, only card punching will occur. 

In the event that any errors occur during computing they will be noted on 
the typewriter. The output data can be edited by dropping in punched 
title cards at appropriate points in the deck of cards before printing on an 
IBM 407 accounting machine. The data output cards can also be used as 
input cards for other programs or may be stored as such. It is strongly 
recommended that care and attention to detail be given to the preparation 
of the punched cards. The lack of a decimal point or failure to fill in 
unused fields with zeros will lead to difficulties. 

The program has been processed on two separate IBM 1620 systems 
at the Jet Propulsion Laboratory and an independent check was run by 
Computermat, Inc.* using the source program and data presented in this 
paper. 


^Computermat, Inc., 1111 Wilshire Blvd. , Los Angeles 17, California. 
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Table II. GAS LIQUID CHROMATOGRAPHY RETENTION DATA 



FU 

PIG 

P 

T 

Z 


65. 2200 

10. 60000 

729. 0000 

21. 0000 

2. 0000 


SWT 

RTAIR 

D 

TC 

N 


4. 4000 

. 6760 

1. 0700 

191. 1000 

7 


FSTP 

PT 

PI 

PR 

PF 


56. 606 

14. 096 

1277. 176 

1. 751 

2.069 


PGCF 

VP 

TR 




. 709 

18. 646 

21. 542 




RT 

W 

VGP7 

VRO 

VG7 


CRT 

TP 

VGP 

CAY 

VG 

(HEXANOL) 

2.515 
2. 177 

. 265 
1441. 133 

33. 757 
35. 192 

154. 849 
36. 120 

19.860 
20. 704 

(HEPTANOL) 

3. 505 
3. 167 

.315 
1980. 956 

49. 108 
51. 197 

225. 267 
52. 546 

28.891 
30. 120 

(OCTANOL) 

4. 915 
4. 577 

.495 

1577.453 

70. 972 
73. 990 

325. 560 
75. 940 

41. 755 
43. 530 

(NONANOL) 

6. 945 
6. 607 

. 645 
1855. 004 

102.450 
106. 807 

469. 953 
109. 622 

60.274 

62.837 

(DECANOL) 

9. 765 
9.427 

.950 
1690. 508 

146. 179 
152. 395 

670. 538 
156. 411 

86. 000 
89. 657 

(UNDECANOL) 

13. 895 
13. 557 

1.270 

1915.268 

210. 220 
219. 159 

964. 303 
224. 935 

123. 677 
128. 937 

(DODECANOL) 

19. 735 
19. 397 

1. 650 
2288.897 

300. 777 
313. 568 

1379. 700 
321. 832 

176. 954 
184. 479 
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The conditions for the above separation were: the column was 6 ft in 
length with 25% wt/wt of Carbowax 20M on a chromosorb P 30/60 mesh 
solid support. The total weight of liquid phase was 4.40g. The helium 
carrier gas flow of 65. 22 ml/min was measured at 21. 0°C and 729. 0 
mm. The column temperature and inlet pressure were 191. 1 C and 
10.60 lb/in. ^ gauge, respectively. The recorder chart speed was 2.00 
cm/min, and the density of the liquid phase at 191. 1°C is 1. 070 g/ml 
(Ref. 6). The proper organization of the above measurements is shown 

below. 





24 

30 

33 

42 

48 

54 

j 60 


10-60 

7 2 ,9 , • ,0, 

2 , l , ■ ,0, , 

2, -.0,0, , 

4, -,4.0. , 

0. ■ .6,7,6, 

1 , • ,0, 7,0, 

1 ,9. 1 , - , I , 

0 , 7 1 , , , 


0, • ,5,3, , 

7 , • .0, 1 , . 

0, ■ ,6 ,3. . 

9. • .8 , 3 , . 

0. • .9 ,9, , 

i . 3 . • , 0 , 9 , 

1 , ■ , 2 ,9 . , 

. , , , , 

. 1 , . . 

1 ,9. • ,5,3 , 

1 , ■ .9,0, , 

2 , 7, ■ , 7 ,9, 

2. -.5.4, , 

3.9, ■ ,4,7 , 

3 , • .3,0, , 

0 • .0 . . , 

0. • .0, , . 

- 

. 1 ■ ■ ^ 


Data for test case 
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